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Using Riemann-Stieltjes methods for integrators of bounded p- variation we define a pathwise in- 
tegral driven by a fractional Levy process (FLP). To explicitly solve general fractional stochastic 
differential equations (SDEs) we introduce an Ornstein-Uhlenbeck model by a stochastic inte- 
gral representation, where the driving stochastic process is an FLP. To achieve the convergence 
of improper integrals, the long-time behavior of FLPs is derived. This is sufficient to define 
the fractional Levy-Ornstein-Uhlenbeck process (FLOUP) pathwise as an improper Riemann- 
Stieltjes integral. We show further that the FLOUP is the unique stationary solution of the 
corresponding Langevin equation. Furthermore, we calculate the autocovariance function and 
prove that its increments exhibit long-range dependence. Exploiting the Langevin equation, we 
consider SDEs driven by FLPs of bounded p-variation for p < 2 and construct solutions us- 
ing the corresponding FLOUP. Finally, we consider examples of such SDEs, including various 
state space transforms of the FLOUP and also fractional Levy-driven Cox-Ingersoll-Ross (CIR) 
models. 
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1. Introduction 

In this paper we consider (stationary) solutions to SDEs of the form 

dX t = p,{X t ) dt + a(X t ) dLf , teR, (1.1) 

where L d is a fractional Levy process (FLP) of bounded p-variation for p < 2 and [x and 
o~ are appropriate coefficient functions. Applying pathwise Riemann-Stieltjes integration 
for functions of bounded p-variation, we solve such equations by constructing explicit so- 
lutions. The basic model will be a fractional Levy-Ornstein-Uhlenbeck process (FLOUP) 
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introduced by the stochastic integral representation 

L d t ' x = f e- x[t - u) AL d ul teR. 

J — oo 

We further show that this is the unique stationary pathwisc solution of the corresponding 
Langevin equation 

dCf x = -\Cf x dt + dLf, teR. (1.2) 

Using this relation we will consider SDEs of the form (1.1) and impose assumptions on 
the coefficient functions /i and a, under which solutions can be constructed by monotone 
transformation of £ d ' x . 

Although our paper is purely theoretical, we are aiming at applications where positive 
solutions of (1.1) are of interest. An approach, developed in [1] for SDEs driven by FBM, 
can be modified to SDEs driven by FLPs. On the other hand, a squared FLOUP is 
positive and a solution to the SDE 

dX t = -2\X t dt + 2^/\xT\dLf, teR. 

We will discuss various examples with different state spaces and different \x and a. We will 
also present some properties of the respective solutions, also concerning the stationary 
distribution. 

Our paper is organized as follows. Section 2 considers FLPs and pathwise integration. 
Section 3 introduces the FLOUP as a pathwisc improper Ricmann-Sticltjes integral and 
shows that it is the unique stationary pathwise solution of the corresponding Langevin 
equation. Moreover, we calculate its autocovariancc function and show that the incre- 
ments of an FLOUP exhibit long-range dependence. Section 4 mainly extends Buchmann 
and Kliippelberg [1] from fractional Brownian motion to FLPs and states structural con- 
ditions for the coefficient functions \i and a, which guarantee an existence (and unique- 
ness) theorem. Section 5 provides examples and simulations. The Appendix reviews the 
Riemann-Stieltjes analysis via p-variation. 

The following notation will be used throughout. We always assume a complete proba- 
bility space (J7,.F, P). We denote the ^-measurable real functions by L°(fl), the Hilbert 
space of square intcgrable random variables by L 2 (£l), the vector space of continuous real 
functions on A CM by C°(A) and by || • |^ p the supremum norm. Furthermore, Lip(A) 
and C (A) are the spaces of real functions on A, which arc Lipschitz continuous on com- 
pacts and continuously differentiable, respectively. The spaces of integrable and square 
intcgrable real functions are denoted by L 1 (R) and L 2 (R), respectively. When speaking 
of a two-sided Levy process we mean the following: given two independent copies of the 
same Levy process, L 1 and L 2 , we take 

Lt ■= £(l{t>o} + £_ t _l{t<o}) teR. (1.3) 

The Dirac measure in 1 we denote by 5i . Finally, for — oo < 6 < a < oo we set [a,b] := [b, a] . 

Integrals throughout this paper are considered in the Riemann-Stieltjes sense, if not 
stated otherwise. 
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2. Fractional Levy processes and pathwise integration 

Fractional Levy processes (FLPs) were introduced as a natural generalization of the 
integral representation of fractional Brownian motion (FBM). We shortly review the main 
properties of FLPs, see [12], Section 3, for details and more background. For notational 
convenience we work with the fractional integration parameter d G (-|,|) instead of 
the Hurst index H, where d = H — i. Because we are only interested in long memory 
models, we restrict ourselves to d <G (0, A). Furthermore, we only consider FLPs with 
existing second moments. Analogously to [11] for FBM we choose (like Marquardt [12]) 
the following definition. 

Definition 2.1. Let L = (X f ) tS R be a zero-mean two-sided Levy process with E[L(1) 2 } < 
oo and without a Brownian component. For d £ (0, ^ ) we define 

Lt-Y^^J [(t-a)i-(-s)i]L{ds), ten. (2.1) 

We call L d = (Lf)teR a fractional Levy process (FLP) and L the driving Levy process 
ofL d . 

The integrals above exist in the L 2 (f2)-sense; see [12], Theorem 3.5, for details. 

Recall that, by the Levy-Ito decomposition, every Levy process can be represented as 
the sum of a Brownian motion and an independent jump process. The Brownian motion 
gives rise to an FBM, which has been studied extensively; see, for instance, [13] for general 
background or [1] in the context of the present paper. 

The next result ensures that there is, in fact, a modification of (2.1) that equals a 
pathwise improper Riemann integral and gives first properties. 

Proposition 2.2 ([12], Theorems 3.7, 4.1 and 4.4). Let L d be an FLP with de 
(0, ^). Then the following assertions hold: 

(i) L d has a modification that equals the improper Riemann integral 

[(t~s) d - 1 -(-s) d - 1 }L(s)ds, teR. (2.2) 



T(d) 



Furthermore, (2.2) is continuous in t. 

(ii) For s,tgR we have 

(iii) L d has stationary increments and is symmetric, i.e., (Ll t )t £ R = (— Zj)( 6 r. 
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From now on, we always work with the modification of Proposition 2.2(i). 

Next we define integration with respect to FLPs. As has been shown in [12], Theo- 
rem 4.10, FLPs may not be semimartingales, and integration in the L 2 (f2)-sense has been 
developed in [12], Section 5. Theorem 4.4 in [12] also shows that FLPs are only Holdcr 
continuous up to the fractional integration parameter d and not to the Hurst index H 
as in the case for FBM. Therefore, pathwisc Riemann-Sticltjes integration by Holder 
continuity does not make sense for SDEs. On the other hand, using an approach like 
Young [17] based on p- variation of the sample paths, integration in a pathwise Riemann- 
Stieltjes sense can be defined; for details see the Appendix. This means we have a chain 
rule and a density formula provided the integrator is of bounded p- variation for p £ [1,2). 

We recall the definition of p- variation over a compact interval [a, b] C R. Let / : [a, b] n- 
R. We define for < p < oo the p- variation of / as 

n 

v p (f, [a,b]) ■= sup^l/(^) " f(xi-i)\ p , (2.4) 

K 

1 — 1 

where the supremum is taken over all subdivisions re of [a, b]. If v p (f, [a, b]) < oo, then we 
say that / is of bounded p -variation on [a, b]. We will further call an FLP L d of bounded 
p- variation if it is a.s. of bounded p-variation on compacts. 

Let L d be an FLP of bounded p- variation, d £ (0, |) and p £ [1,2). For AcRwc define 

W C ™(A) := {/ G C°(A) : v p {f, [s, t}) < oo V[s, t] C R}. (2.5) 

Then we define for every stochastic process with sample paths H £ 2U^ on (M) a.s. and for 
p, q > with + q^ 1 > 1 the integral 

b 

H s dL d s: -oo<a<fc<oo, (2.6) 

pathwise in the Riemann-Stieltjes sense. 

As stated in the Appendix the integral in (2.6) always exists on finite intervals [a, b]. 
We consider also improper integrals, where a = oo or b = — oo. The existence of the tail 
integral has then to be justified. 

For example, FLPs, where the driving Levy process is of finite activity, are of bounded 
p- variation for all p> 1; cf. Theorem 2.25 of [12]. 

3. Fractional Levy— Ornstein-Uhlenbeck processes 

We introduce fractional Levy-Ornstein-Uhlenbeck processes (FLOUPs) as improper 
Riemann-Stieltjes integrals and prove that they are stationary solutions of the Langevin 
equation (1.2). To show the existence of the improper Riemann-Stieltjes integral, we first 
need some knowledge about the long-time behaviour of FLPs. A similar result considering 
t — > oo has been proven by Muncya Matsui (personal communication) . 
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Theorem 3.1. Let L d be an FLP, d £ (0, ^). Then for all a> d+ \ we have 

\L d t \ 



lim 

t— > — OC 



a.s. 



(3.1) 



Proof. Without loss of generality we can assume that t < 0. By the law of the iterated 
logarithm (LIL) for Levy processes (cf. [14], Proposition 48.9) we find a random variable 
T and a constant M > such that a.s. for all s < T 



|L( S )|<M(2| S |loglog| S |) 1 / 2 . 



(3.2) 



We can always make T smaller and so we choose T < — e. For any such path we can 
assume that t <T and estimate 



1 



1 1 



< 



\A a r(d) 
i i 



— oo 

t 



r(d) \t\ 

Therefore, it suffices to show that a.s. 



[(t-*)t "(-«)+ ]£(*)d* 
(f [(- S ) d - 1 -(t-s) d - 1 }\L( S )\d S + /'W-^OOIds 

\J~oo Jt 



(3.3) 



lim — — 

t^-oo \t\ a 



[(-sr-i-(t- s r- i ]\L(s)\ds=o 



and 



^jT(-.)^(-)id-=o. 



(3.4) 



(3.5) 



We start with (3.4). Using the LIL we get an upper bound of (3.3) as follows: 
^ f l(-a) d - 1 -{t-8) d - i ]\L(s)\da 



< 



\t\° 

M\t\ 
e\t\« 



- (-|t| - ,s)^ 1 ](2| S |loglog| S |) 1 /2 ds 

[(-e- 1 |t|«)««- 1 -(-|t|-e- 1 |i|«)«'- 1 ] 
x (2c- 1 |t||u| log log(e- 1 |t||w|)) 1 / 2 du, 



(3.6) 



where we have used in the last line the change of variable e 1 \t\u = s. Now note that for 
large |i| and \u\ > c 



|u| loglog(e- 1 |*|| W |) = \t\\u\ log(log(e- 1 |*|) + log 
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< \t\\u\ loglog |«| + |.||u| log(l + log |u|). 

Combining (3.7) with \a + b] 1 / 2 < \a\ x / 2 + lb] 1 / 2 for a, b £ R we get an upper bound for 
(3.6) by 



M{2c- 1 \t | loglog \t\f/ 2 
e\t\ a ~ d 

M(2e~ 1 \t\) 1 / 2 



e|t|" 



-d 



[{-c^uY- 1 - (-1 - e- 1 «) ,, - 1 ]|«| 1 / 2 du 

oo 

[(-<T l u) d - 1 - (-1 - e- 1 u) tf - 1 ](|tx| log(l + log \u\)) 1/2 du 



(3.7) 

+ jT^"" 1 "^" 1 - (-i+c- v-^ioga+iogM)) 1 ^. 

By a binomial expansion we get (e _1 it — l) d_1 = (c~ 1 u) d ~ 1 — (d — l)(e _1 ii) d ~ 2 + 0(u d ~ 3 ) 
and, therefore (writing a(u) ~ .(it) for lim.u_i.oo a{u)/b{u) = 1), 

[(e" 1 ^)^ 1 - (-1 + c- 1 w ) d - 1 ]( w log(l + log M)) 1/2 

(3.8) 

-(d-^Ce-^^-V-^CloglogW) 1 / 2 , 

which ensures the existence of the two integrals in (3.7). 
Letting t — > — oo, we obtain (3.4). Next we calculate 

1 I ' \d-\\T I „\l A„ 1 / I „\d-\\ T ( „M A „ / / „\_ — 1| 



^ Jf (-_)-* |L(.)|d_ = ^ jf (-_)-* |L( s )|da+|^ y^ (-.)- i |X(.)|d*. 
The second term tends to zero as t — > — oo, and we consider the first: 

— J (- s )^|L( s )|d.<— J (- S ) rf - 1 (2| S |loglog| S |) 1 /2 ds 
< M(2\t\ loglog I*!) 1 / 2 [ T , )d -i ds 
_ i\/(21oglog|t|) 1 /2 |T| d M(2 loglog Itl) 1 / 2 

d|t|«-(«*+l/2) _|i|«-l/2 ' 

Letting £ — ^ — oo, we get (3.5) and therefore the assertion. □ 



Theorem 3.1 ensures the existence of the improper Riemann-Stieltjes integral. 
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Lemma 3.2. Let L d be an FLP, d G (0, |) and A > 0. Then for — oo < a < oo 

f e Xs dL d , t>a, (3.9) 

J a 

exists a.s. as a Riemann-Stieltjes integral and is equal to 

e xt L d -c Xa L d -\ [ L d e Xs ds. (3.10) 



Furthermore, the function (a, oo) — > R defined by 1 1— > J*e Xs dL d is continuous. 

Proof. From Theorem 3.1 - 
that for ui G SI \ N we have 



Proof. From Theorem 3.1 we know that for all a > d + \ there is a null set 2V C SI such 



l im ^1=0 (3.11) 

t^-oc. \t\ a ^ ' 

and, hence, for all u> G \ JV and < > a, the integral X^(w)e A " dit exists as a Riemann- 
Sticltjcs integral. For a compact interval [a,t] this is clear. Now consider a~ —oo. It 
suffices to show that J_ L d (uj)e Xu du exists for T < — 1. This follows from the inequality 

T rT \y d/ \\ r-T 



L d {uj)e Xu du 



|L "^ )l e A "k| a d M <C / e Au |u| Q du 



for some constant C > 0, and the integral on the right-hand side exists for R — oo. 
Similarly, 



lim c Xa L d a {uj)=Q. (3.12) 

a— >■— oo 



Now it follows by Wheeden and Zygmund [16], Theorem 2.21, that (3.9) also exists as a 
Riemann-Stieltjes integral and is equal to (3.10). Since (3.10) is continuous in t for all 
t> a, the result is proven. □ 

Now we are ready to define the central object of this paper. Recall that all integrals 
are Riemann-Stieltjes integrals based on Lemma 3.2. 

Definition 3.3. Let L d be an FLP, de (0, ±) and A > 0. Then 

Cf X := f e- x{t - s) dL d s , t&R, (3.13) 

J — oo 

is called an FLOUP. 



Before returning to the Langevin equation in connection with the FLOUP, we present 
some distributional properties of £ d ' X . With a little effort one can prove that C d > x is 
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stationary, i.e., for all t\ < • • • < t m , m£N, h&M., 

(£f; A ,...,0 (3-14) 

For details see [5], Lemma 6.1.3. 

Although we mainly concentrate on Riemann-Stieltjes integrals, there exist several 
results based on integrals in the L 2 (fi)-sense that we can use. Fractional integration 
can be considered as a transformation of classical Ricmann-Liouvillc fractional integrals, 
which are defined for < a < 1 by 

(I™f)(x) = — f{t){t-x) a - x dt and 



Via) J„ 

{i+s){x) = yj-r / fmt-xr-'dt, 

r (") J-oo 

if the integrals exist for almost all ieR. This is, for instance, the case if / € L' P (M.) with 
1 < P < h- The following result is a Riemann-Stieltjes version of Theorem 3.5 of [12]. 

Proposition 3.4. Let C d ' X be an FLOUP driven by an FLP L d of bounded p-variation, 
d £ (0, i), A > and p > 0. Then its finite- dimensional distributions have a characteristic 
function 



E 



expjf>,4 A j =exp|^^E^( J - e " Afe "' )l ^>-})( s )J ds }> 



3- 

Ul,...,U m €. 



for — oo <*!<•••< t rn < oo and ^l(u) := J R (e mx — 1 — iux)v(dx), where v is the Levy 
measure of L. Furthermore, for every t £ K, the random variable C d ' X is infinitely divisible 
with a characteristic triple given by {pf^ diX ,Q, v^ d ,x)i where 

7^.* = - / / (^ e " A(t ~ ^t>.})(*)sl { | ( i-e-»c*-oi it >. } ) Wa ,|>i} du(x) ds, (3.16) 
Vcd,x{B)= [ [ l B {{L c lc- x{t -h {t >. } )(s)x)dv(x)ds VBorel sets B in R. (3.17) 



Proof. For simple functions, the Riemann-Stieltjes integral and the L 2 (f2)-integral agree 
a.s. (see [12], Proposition 5.2). Now approximate the function e _A ( i_s 'l{t> s } by simple 
functions. While a.s. and L 2 (f2)-convergence of the Riemann-Stieltjes sums imply both 
convergence in probability, the integrals are equal in probability and thus in distribution. 
Therefore the result follows as in Theorem 3.5 of [12]. □ 

We now turn to the second-order properties of an FLOUP. Chcridito, Kawaguchi and 
Maejima [2] present details concerning the long memory property of an OU process driven 
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by FBM. Similarly, we shall show that the increments of the FLOUP exhibit long-range 
dependence. First, however, we need the following result (see also Proposition 4.4 of [9] 
and Proposition 5.7 of [12]). 

Theorem 3.5. Let L d be an FLP of bounded p-variation, d <G (0, |), A > 0, p > 
and f,g:W—>M. with |/|, \g\ e 2U™ n (IR) for p- 1 + q- 1 > 1, such that J R /(s)dLf and 
/h ff( s ) ^L d sxist as Riemann-Stieltjes integrals. Then we have 



E 



f{t)&L d t / g{s)ALi 
Jr. 

Y(l-2d)E[{L{l)f 
T(d)T(l-d) 



(3.18) 

/(%(«)|t- al 21 *- 1 da <tt. 



Proof. The proof follows again by using approximating simple functions and the fact 
that 

\t - ur-\s - du = r( ^ (1 "^ l* - ^l 2 ^ 1 (3-19) 

r(i - d) 

for i.sgt, which can be found in Gripenberg and Norros [6], page 404. □ 

Now we have everything together to derive the covariance structure of an FLOUP. 
The lengthy calculation works in a manner similar to that of Theorem 2.3 of [2]. The 
asymptotic is up to a multiplicative factor the same as in the FBM case. 

Theorem 3.6. Let L d be an FLP, d <G (0, \), A > and C d ' X the corresponding FLOUP. 
Then for N € No and for fixed t gK we have as s — > oo 

Cov(£f\£^) 



r(l - 2d)E[L{\)] 2 



2d-2N-l\ 



Now it is clear that the increments of an FLOUP exhibit long-range dependence in the 
sense of a non-summability property of the autocovariance function. 
We now return to the Langevin equation presented in (1.2). 

Theorem 3.7. Let L d be an FLP, de (0, f) and A > 0. Then the unique stationary 
pathwise solution of (1.2) is given a.s. by the corresponding FLOUP 

C d ' X = I c- x{t - u) AL d „ teR. 
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Proof. From Lemma 3.2 we know that f e~ A ( t_u ) dL d exists for all t G R a.s. as a 

J — oo w 

Ricmann-Stieltjes integral. We fix s G R and consider the pathwisc SDE 



Cf x = £ s - A I C d u ' x du + L d - L d , s < t, 

J S 



(3.20) 



where £ s := /* e A ( s u ) dL d . Obviously, £ s <G L°(Q). By arguments similar to those in 
the proof of [2], Proposition A.l, we obtain 



[e As f c- x{s - u ^dL d u + f e Xu dL d \= f e~ x{t ~ u) dL d , te 

L J —oo Js ) J —oo 



is the unique pathwise solution of (3.20) and, therefore, by (3.14) a stationary solution 
of (1.2). 

On the other hand, let (X t )t£R be a stationary solution of (1.2). We show that 



(X t ) t £g. — (Cf X )teM holds for almost all u G fi. Set A := {w G : (X t (w)) te a 7^ 



(£f A (w)) teR } and assume that > 0. For we A fix (el with X t (w) ^ Cf ( 

Then we have for s < t by (3.20) 



: t d ' A l 



+ \X t - £ 



At I „As 



c AS X, 



c A " dL d , 



-X(t-v) ^d 



-X(t-u) dL d 



= e~ xt e Xs \X,-C 



d,Xi 



where we supressed the chosen to for simplicity. Since A > and s — > — 00 we conclude 
that \X s (u>) — C d ' X (uj)\ — > 00 for s —> —00. Therefore, on A we have \X t — C d ' X \ —> 00 
for t — > —00. For a given K > we define w-wise the random variable T : A — > R with 
|X t - £?' A | > pf^y for i < T on A. Hence, 

E\X t - Cf x \ > E[\X t - Cf x \ht<T}U + E[\X t - £ d ' x \l {t>T} l A ] > j^P({t <T}n A). 

Furthermore, we know that {t < T} P\ A C {s < T} P\ A for s < t. Choosing a sequence 
(t n )n£N of real numbers with lirrin^oo t n = —00 we get by continuity of P 



lim P({t n <T}nA)=P[ \\{t n <T}r\A) =P(A). 

n— >oo \ — / 

Putting everything together we arrive at 

lim E\X tn -C d t ' X \> lim -£—P({t n <T}nA)=K. 

n— too n— )-oc Jr[J±) 

Hence, lim TWOO E\X tn — C t ' | = oo. However, we have now 

E\X tn I = E\X tn - Cf x - (-Cf x )\ > E\X tn - C d ' x \ - E\C d '\ 



(3.21) 
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where E\C d ' X \ is independent of t n . Thus, lim„_ >00 E\X tn | =00 and, by stationarity, 
E\X t \ = 00 for all t G E. However, we also have for fixed s <t 

lim (X t - C d t ' X ) = lim C - At ( c Xs X s - f c Xu dL d u )=0 a.s. (3.22) 

t-Kx> t^oo \ 

Hence, by stationary X t = C d,X , but E\C d,X \ < 00, which is a contradiction and, thus, we 
conclude that P(A) = 0. □ 

The following Ornstein-Uhlenbeck operator will be used in the next section to obtain 
solutions to SDEs with different starting values. The operator here modifies the starting 
value of the FLOUP and the lemma shows that this modified process still solves the 
Langevin equation. 

Definition 3.8. Let L d be an FLP, d g (0, \), X > and C d x the corresponding FLOUP. 
We define the Ornstein-Uhlenbeck operator by 

£ x {L d ,- r ):RxR — >C°{R), 

(3 23) 

(r, z) 1 y Cf X - e ~ A ('-^> A + C - X ^z. 

It is immediate from this definition that £ A (L d ,r, z) = z a.s. for (r, z) £ M 2 . 
The next lemma shows that L d transformed by the Ornstein-Uhlenbeck operator still 
satisfies the Langevin equation; its proof follows directly by the definition. 

Lemma 3.9. Let L d be an FLP, d g (0, \), A > and C d ' X be the corresponding FLOUP. 
For a continuous process I := (lt)teM ^ e identity l t = £ A (L d , r, l T ) holds for all t, t g R if 
and only if 

dl t = -\l t dt + dL d , teM. (3.24) 



4. Solutions of fractional SDEs by state space 
transforms and proper triples 

In this section we start with SDEs driven by FLPs. Using pathwise integration we must 
solve for almost all w g f2 a deterministic integral equation. Consequently, we build on 
an extensive theory starting with the seminal work by Young [17]. We also recall that for 
Brownian motion the pathwise approach goes back to [4] and [15] leading to the Fisk- 
Stratonovich integral. Required is that [i is Lipschitz-continuous and a g C 2 (R) with 
bounded first and second derivatives. Readable accounts on the history can be found in 
[8] and in Ikeda and [7]. 

Regularity assumptions of sample paths of the driving process like Holder continuity 
or bounded p- variation for p < 2 have been considered by Lyons [10]. We shall work in 
the framework of p- variation, however, to go beyond the work of Lyons, who proves only 
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existence and uniqueness theorems under certain Lipschitz assumptions on the coefficient 
functions and gives no analytical form of the solution. 

The approach by Zahle [19] is indeed comparable to ours, where explicit solutions 
can be given under differentiability and Lipschitz conditions on the coefficient functions. 
Most of her results can be applied to SDEs driven by an FLP of bounded p- variation for 
p < 2. We believe that the contribution of our work is two-fold. First, our assumptions 
are easy to verify and, second, we are able to present analytic solutions to SDEs of the 
form 

dX t = (a\X t \~< + pX t )dt + o-\X t \~< dLf, tel. (4.1) 

In this situation we cannot apply the results of [19], since the volatility coefficient docs 
not match the required differentiability assumption. Lyons [10] provides us at least with 
an existence theorem, but gives no closed form solution. 

Aiming at solutions to similar SDEs, driven however by FBM, Buchmann and 
Kliippclberg [1] presented a theory that can be modified to cover SDEs driven by FLPs. 
The idea is to present solutions to SDEs like, for instance, (4.1) as monotone transfor- 
mations of the FLOUP. The question we shall answer is, given an SDE 

dX t = n{X t ) + a(X t )dLt (4.2) 

for specific /z and a, which monotone transformation of the FLOUP is a solution to (4.2)? 
First we have to establish certain regularity conditions on the coefficients fi and a . 

Definition J^.l. (i) A triple (I, fx, a) is called strongly proper if and only if it satisfies 
the following properties: 

(PI) / = (a, b) C R is an open interval, where — oo < a < b < oo and n, a G C°(I). 
(P2) There exists a strictly decreasing ip, absolutely continuous with respect to the 

Lebesgue measure such that ip — fi/a on I\ Z(o~) where Z(o~) are the zeros of a, 

and 

lim ip(x) = — lim ip( x ) = — °°- 

(P3) There exists A > such that aip' = A holds on I Lebesgue-a.e. 

(P4) The inverse function ip -1 :R — > tp~ 1 (M) = I is differ entiable and (^ _1 )' £ Lip(R). 

(ii) We call the triple (/, /x,<t) proper if only (P1)-(P3) are satisfied. 

(iii) The interval I is called the state space, the unique constant A > in (P3) is called 
the friction coefficient (FC) and the unique function f :K — > I = f(R), f(x) := ^~ 1 (— Ax), 
is called the state space transform (SST) for (I,fj,,a). 

Condition (P4) differs from the H -proper assumption required in [1] , because we work 
with p- variation instead of Holder continuity. 

As pointed out in [1], > tp(I) = M is by (P2) strictly decreasing and a.c. differ- 

entiable on / with ip' < 0. Condition (P3) implies that Z(a) and Z(ip') have Lebesgue 
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measure zero. Also we have that a is non- negative and 1/tr £ Cc(I), where Cc(I) denotes 
the locally integrable functions on 7; 7\ Z(a) is dense and open in 7 by (PI). Therefore, 
the equality pL = ai\> extends to 7. It follows that ip and A are uniquely determined by /z 
and a. 

As can be seen from Definition 4.1 our coefficient functions are only defined on the 
interval 7, which can be any interval in R. To account for this situation we need to specify 
what will be understood to be a solution to an SDE. 

Definition 1^.2. Let L d be an FLP of bounded p-variation, p£ [1,2) and d £ (0, |). 
Suppose that ICR is a non-empty interval and n,cr £ C°(I). We refer to a stochastic 
process X := (Xt)teM as a pathwise solution of the SDE 

dX t = /i(X t )dt + o-(X t )dL d , t£R, (4.3) 

if for almost all sample paths the following conditions are satisfied: X £ 2U™ n (R) and the 
image of X is a subset of I such that for s <t: 

(51) u o X is a.s. Riemann-Stieltjes integrable with respect to L d on [s,t]; 

(52) The following integral equation holds: 

X t -X s = j n(X u )du + J^ o-{X u )dL d u . 

The space of all solutions of (4.3) is denoted by S(I,fj,,o~,L ). 

We consider now an SDE as given in (4.3). If we assume that the triple (J, /i,cr) is 
strongly proper with SST / and FC A, we define the following operator 

A^ A (L d ,.,.):Mx/^C°(R), (t, z) i > f(£,^(L d , t, f~ 1 (z))) (4.4) 

with Ornstein-Uhlenbeck operator as in Definition 3.3. We also remark that 

X{'\L d ,r,f(C^)) = f(Ct' x ), t£R. (4.5) 

Before we state our main results we state the following technical lemma, which will be 
needed in the proofs. 

Lemma 4.3 (Version of Lemma 3.2 [1]). Let (I, H,cr) be a strongly proper triple 
with the corresponding SST f. Then f £ C 1 (R) with derivative f = a o f. Also f £ 
C 1 (/\Z(cr)) with (/ -1 )'(a;) = l/cr(x) for all x £ L \ Z(a) . 

Next we state the existence theorem. Let M(f2,7) denote all mappings from f2 into I. 

Theorem 4.4. Let L d be an FLP of bounded p-variation, p £ [1,2) and d £ (0, |). If 
(I,fi,a) is strongly proper with SST f and FC A > 0, then 

{X f ' x (L d , T , W) : t £ R, W £ M (fi, 7)} C S(I, n, a, L d ). 
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Proof. Because we consider pathwise solutions we can w.l.o.g. assume that W = z a.s. 
for some z £ I. Now fix r € R and z E I and define 

U'.= ${L A ,T,f-\z)) and Y t :=Xi' X {L d ,T,z), teR. 

We show that Y = (It)tgH € S(I,/j,,a,L d ). Obviously V takes only values in /. Since 
/ e C*(R) and I is of bounded p- variation, we know that Y = fole 2B™ n (R). With the 
chain rule from Theorem A. 2 we get for s,t 6 R 

Y t -Y s = f(l t )-f(l s )= f f'(l u )dl u , s<t, (4.6) 



since I solves (3.24), 

l u = l s -X I l v dv + Lt - L d s , s<u. (4.7) 



The Riemann— Stieltjes integral is additive with respect to a sum of integrators, if the 
Ricmann-Stieltjcs integrals exist separately for each integrator. This is true in our case, 
because J l v dv is of finite variation and L d is of bounded p-variation. Since f'(l u ) is 
continuous and also of bounded p- variation, (4.6) and (4.7) imply for s,teR 

Y t -Y s = jj'{l u )d(-\j^ l v dv^j+ Jj'{l u )dL d ul s<t. 



Furthermore, J u l v dv is differentiable and f'(l u )lu is continuous as a function of u and, 
thus, we get by the density formula for Riemann-Stieltjes integrals on compacts for 
s, t e R 

Y t -Y. = -X [ f'(l u )lu du + f f'(l u ) dL d u , s < t. 

J S J S 

From Lemma 4.3 we obtain /' = a o f, hence a o Y = f o I e W c ° n {R). By Defini- 
tion 4.1(P3) and the interpretation following this definition we find that <r/ _1 = —aip/\ — 
— /i/A. This yields 

Y t -Y s = -\J^ f'{l u )l u du + j f{l u )dL d u 

= -\j a(f(l u ))l u du + J a(f(l u ))dL d u 

! ft 

^{Y u )du+ / a(Y u )dL d u , s<t, 



where we used in the last line the equality <r(f(l u ))l u — a(Y u )f 1 (Y U ) = —fj,(Y u )/X. Fi- 
nally, we have Y € S(L, /i, a, L d ) . □ 
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The following result ensures uniqueness under natural conditions. 

Theorem 4.5. Let L d be an FLP of bounded p- variation, p£ [1,2) and d£ (0, |). Let 
also (I,/j,,a) be strongly proper with SST f and FC A > 0. Furthermore, assume that 
Z(cr) = 0. Then 

{X f ' X (L d , t,W):t£R,W £ M(Q, 7)} = S{I, n, a, L d ). 

Proof. From Z{a) = 0, we know by Lemma 4.3 that / £ C^R) and = l/a(x) 

for all x £ L. Let X £ S(L,n,a,L d ). From Definition 4.1 we know that X £ 2B£ on (R) a.s. 
and from (Z" 1 )' £ Lip(7) we get by the chain rule from Theorem A. 2 for s,t £ R 

r i {x t )-r i {x s )= f r 1 {x u )dx u = f-i- dx u , 8 <t. (4.8) 

Js Js 0\A u J 

Since X € 5(7, /i, cr, L d ), we know that for s, u £ R 

/i(X„)d« + y cr^dL*, s<w. (4.9) 

Now J u fi(X v ) dv is of finite variation and by the density formula of Theorem A. 3 the 
integral J" cr(X v ) dL d is of bounded p-variation as a function of u. By putting (4.8) and 
(4.9) together and using again Theorem A. 3, we get for s,tgR 

f-\x t )-r\x s ) 

t -j^ f f U \ 1 



i (/"" (x - )d ") + /^ d (/ %(x - )di 



du + L?-Zj, s<f, 



since (I, /i,cr) is proper, V'( a; ) = M^X* 7 ^)) 1 an d VK^) = — A/ (a;) hold for all x £ I. 
Thus, 



f~\X t ) - f-\X s ) = [ ^Mdu + Lf-Lt 

Js CT(^«) 

^(X„)du + L t d -Lf 



-Ay /"^X^du + if-Lf, s<i. 
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Hence, f~ 1 (X) is a solution of (3.24). Fixing relwe see by Lemma 3.9 that f^ 1 (X) = 
£, x {L d ,T,f- 1 (X T )) and, finally, X = Xf< x (L d ,T,X T ). □ 

The next corollary covers the important case of a stationary solution. 

Corollary 4.6. Let L d be an FLP of bounded p-variation, p£ [1,2) and d£ (0, ^) and 

C d ' X be the corresponding FLOUP. Furthermore, let (I, /i, a) be strongly proper with SST 
f and FC A > 0. Set X t = f(C d ' X ) for t £ R. Then the following assertions hold: 

(i) X is a stationary pathwise solution of the SDE 

dX t =n(X t )dt + a(X t )dL d , teR. (4.10) 

(ii) // Z(a) = 0, then X is the unique stationary pathwise solution of (^.10). 

Proof, (i) From Theorem 4.4 we know that X = X^ x {L d , 0, / (£q ,A )) = f(£ d ' X ) is a path- 
wise solution of (4.10). Furthermore, X is stationary as a transformation of a stationary 
FLOUP. 

(ii) Given a pathwise solution of (4.10), Theorem 4.5 supplies us with a W € M(f2, 1) 
such that Y t = f(c~ Xf W + Cf' X ). If we want Y to be strictly stationary, we must have 
W = a.s. and get Y = X. □ 



5. Examples of SDEs driven by FLPs 

5.1. Examples by means of strongly proper triples 

This section is dedicated to examples, which we illustrate by simulations. For those we 
consider as driving Levy process a compensated Poisson process L e with intensity 9 > 0; 
that is, 

L\ := Pf -~t9, t G M, 

where P is a Levy process with drift 7 = and Levy measure v{dx) = 95i(dx) without 
Brownian component. Of course, we consider this process to be defined on the whole of 
K using (1.3). 

In a first step we simulate sample paths of L e and compute the corresponding FLP L d 
by a Riemann-Stieltjes approximation; that is, we approximate 




500 H. Fink and C. Kliippelberg 

From Theorem 2.55 of [12] we know that the quality of this approximation is 

O^ 1 ' 2 ) + Oin- 1 ' 2 ) + 0(n( 1+2d - 2d2 >/( 2d - 3 >). 

Furthermore, the Poisson-FLP is of finite variation by Theorem 2.25 of [12]. 
Now we use a version of the explicit Eulcr method for the SDE (1.2) 

d£f x = -\£f x dt + dLf, teR, 

to compute sample paths of the FLOUP. We want to remark that all these computations 
are pathwise. Probability comes in only through the underlying paths of the driving Levy 
processes. 

Next we study some examples of solutions to the SDE (4.3) given by strongly proper 
triples. We will mainly draw from structural results of [1] taking into account that their 
-ff-propcr condition has to be replaced by our assumption (P4) in Definition 4.1. 

For the rest of this section, let L d be an FLP of bounded p-variation, p S [1,2) and 

rfe(o,i). 

Example 5.1. As a first example, we consider for parameters a, (3 £ R and a > an 
SDE of the form 

dX t = (a\X t p + f3X t )dt + a\X t p dL?, teR. 

We analyse this SDE by taking the volatility coefficient a : R — > [0, oo) defined by <j(x) := 
<Jo\x\ J for do > and 7 € R as given. The question is now, what drift functions /i and 
intervals I lead to strongly proper triples (I,fi,a) as defined in Definition 4.1? More 
precisely, we want to find elements in the set 

JCi := {(A,/i) e R + x C°(J) : (I,/i,cr) is proper with FC A}. 

Using Proposition 5.5 of [1] we see that only 7 £ [0,1] leads to a non-empty K} a . We 
consider the cases 7 = 0, 7=1 and 76 (0,1) separately. 

Take first 7 = 0. For a triple (I,/j,,a) to be proper we must have that crip' = —A with 
t/j = fi/a. This results in dX t = (a + (3X t ) dt + adLf with state space / = R and SST is 
affinc, more precisely, f{x) = a + (3x for 

If 7= 1, by Proposition 5.6 of [1], the state space can only be either I — (— 00, 0) or 
/ = (0, +00). For / = (0, 00) we get 

£(0,00) = {(|£| )At ) e R+ xC°(I):n(x)=ax + /3xlogx,a€R,j3<0,x£ (0,oo)} 

and the state space transform is f(x) = expjaox — ^} for x € (0, 00). Simple calculation 

ensures that condition (P4) of Definition 4.1 is satisfied, and every element of /cS°'°°' > 
leads to a strongly proper triple. An example of an SDE of this kind for a = can be 
found in (5.6). The case / = (— 00, 0) can be treated analogously. 
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Finally, we consider 7 <G (0, 1). Proposition 5.8 of [1] shows that the only possible state 
space is the whole real line R and 

K* = {((1 - 7)101, A«) e IR+ x C°(I) : n(x) = a\x\~< + 0x, a e R, /3 < 0, x e M}. 



Furthermore, the SST is given by 



f(x) = sign (1 -7)002;- 



(l--y)a x- 



1/(1-7) 



The derivative of / can easily be calculated yielding that only 7 £ [i, 1) leads to strongly 
proper triples. An example for such an SDE (with a = 0) is a fractional Cox-Ingersoll- 
Ross-type model, which is investigated in detail in Section 5.2 (cf. the SDE in (5.3)). 



Example 5.2. We consider the following SDEs with affinc drift. 

dX t = (a + pX t )dt + a(X t )dLf, teR, 

that is, /z:R — > R is defined by fi(x) := a + (3x for a,/)£l. To find suitable volatility 
coefficients and state spaces, wc consider the set 

A£ :={XeR + :3cr eC°{I) with (I,fi,a) is proper with FC A} 

and, if there is a A € A^, we investigate 

K{ l . x = W G C°(J) : (I, ^, cr) is proper with FC A}. 

Proposition 5.1 of [1] implies that there exist I, a with (J, /x,cr) being proper if and only 
if /? < 0. In this case it also follows that I = R and A^ = (0, |/3|]. A FC A = |/3| leads again 
to an affine model, namely <r(x) = a^x for some <7o > 0. 

If we choose an FC A = (1 — 8)\f3\ £ (0, \f3\) for some 5 £ (0, 1), then by Proposition 5.3 
of [1], every a £ "H* n_g\\p\ is of the form 

cr(x) =cn|a + /3a;| A l {;E <_ Q//J} + a 2 \a + Px\ s l {x >_ a/0} 

for some <ti,<t 2 > 0. Setting ft := \p\ 5 /^~ s ) g] /{1 ~ 5) (1 - S) 1 ^ 1 ^ for i = 1,2 the SST takes 
the form 

/(«) = (j - 1{.<o } + (j + / 2 |-| 1/(1 - 5) ) 1 { *> 0} . 

Calculating the derivative of /, we see that a possible proper triple is strongly proper 
if and only if S <G [|, 1). An example of such an SDE is for parameters a £ R and j3 < 
given by 

dA t = (a + /3A t ) dt + a y/\a + pX t \dL$, teR. (5.1) 
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Example 5.3. Consider the SDE 

dX t = -a 1 sm(a 2 X t )cos(a 2 X t )dt-8m 2 (a2Xt)dL d , tel. (5.2) 

This example provides a bounded state space model. Define the triple (I,fi,a) by / := 
(0, ^-), fJ-(x) := — (Ti sin(cr2^) cos(a2^) and u(x) := — sin 2 ((i2a;) where cr\,a 2 > 0. It can 
be shown that this triple is in fact strongly proper with FC A = u\<Ji - More precisely, we 
have tp(x) = o\ cot((722;) and, therefore, 

X t ■= — arccot(-f72£f CTl<r2 ), teR, 

is the unique stationary solution of the SDE (5.2). 

5.2. Fractional Cox— Ingersoll— Ross models 

Whenever positive phenomena are modeled - for instance, interest rates, volatilities or 
default rates in finance - the Cox-Ingersoll-Ross (CIR) [3] model is the most prominent 
model. It is the solution to 

dX t = (a-jXt) dt + a^fx~ t dB t , A o =xo>0, 

where B = (-Bt)te[o,oo) denotes standard Brownian motion, 0,7 € R and a > 0. General 
existence and uniqueness theorems of Brownian SDEs cannot be applied here, because 
the square root is clearly not Lipschitz continuous. However, Ikeda and Watanabe [7], 
page 221, showed that for any X$ = x > there exists a unique non-negative solution. 
We shall consider analogous SDEs driven by FLPs. 

Within the framework of strongly proper triples, Examples 5.1 and 5.2 show that our 
theory only covers CIR models with mean reversion to a = 0. Consider for a, 7 > a 
solution to the pathwise SDE 

dX t = -~fX t dt + aiy\x7\dLd, teR. (5.3) 

Define a(x) '.— alx] 1 ^ 2 , choose Jl(x) = —72; and take / = R. Example 5.1 implies that 
(/, /x,cf) is strongly proper with SST 

^2 

f(x) =sign(a;)— x 2 , 

and, by Theorem 4.4, a stationary solution of (5.3) is given by (/(£f' A ))t 6 R with A = 7/2, 
cf. Figure 1. Obviously, this CIR model takes also negative values. 

A natural non-negative transformation of the FLOUP is given by Z t := (<r£^' X ) 2 for 
t £ R (cf. Figure 2), and, using the chain rule from Theorem A. 2 and the existence of all 
appearing Riemann-Stieltjes integrals, we get 

dZ t = 2a 2 Cf X dC d t X = -2\a 2 (£f> x ) 2 dt + 2a 2 Cf' X dL d t 

= -2XZ t dt + 2a^fz~ t dL d t , teR. 
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t t 

Figure 1. Sample paths of a solution of the Cox-Ingersoll-Ross model (5.3) with Xo = for 
varying a, fixed A = 2.5 and d = 0.35, using two different FLP sample paths, left 9 — 0.5, right 
61 = 2.5. 

Defining now k(z) := —2Xz and i{z) := 2a^/z we have Z as a solution to 

dZ t = n{Z t )dt+i{Z t )dL d , teE.. 

However, the triple ((0, oo), k, l) is not strongly proper, because assumption (P2) of Def- 
inition 4.1 is violated. 

We can now formulate the following general result. 

Proposition 5.4. Let Z x / 2 {L d 1 •, •) be the Ornstein-Uhlenbeck operator from Defini- 
tion 3.8. Then for r £ M and z > i/ie process 

X^- z :=^ /2 {L d ,T,z)j , t G R, 

solves the SDEs 

dX t = -XX t dt + a/^dL t d (5.4) 

and 

dX t = -AXtdt + aV^dif, tel. (5.5) 
In fact, any solution to (5.5) afeo solves (5.4). 

This result is not surprising, because Theorem 4.5 does not hold for the SDE (5.4). 
However, the reverse is not true: a solution of (5.4) does not necessarily solve (5.5), 
because it can be negative. Also the constant process, X t := 0, t E K, solves both (5.5) 
and (5.4). 

Considering a squared Ornstein-Uhlenbeck process leads in the case of a driving Brow- 
nian motion to a CIR model with mean reversion to positive values. This approach does 
not work for pathwise integrals, neither in the FLP case nor for FBM, since the ltd term 
in the chain rule vanishes by finite p- variation for some p < 2 . 
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Figure 2. Sample paths of squared FLOUPs for varying a, fixed A = 2.5 and d = 0.35, using 
the same sample paths as in Figure 1: left 6 = 0.5, right 9 = 2.5. 




A positive process based on Theorem 4.4 is given as a solution to 

dY t = -Xy/Y t \og(Yt)dt + a\Y t \dLf, t e R, Y = y (5.6) 

for A,(T > 0, cf. Figure 3. Example 5.1 states that the triple for state space /= (0, oo) is 
strongly proper and the SST can be calculated as f(x) = c ax . 



Appendix: Riemann— Stieltjes integration 

As mentioned in the Introduction, all integrals in this paper are considered Riemann- 
Stieltjes integrals, if not stated otherwise. That is, for functions /, h : [a, b] >->• R, we take 
the limit of 



n 

S{f,g,K,p) :=^2f{yi)[h{xi) - h(x l - 1 )}, 

i=l 



(A.l) 
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where k = (a;i)i=o,...,n is a partition and p = (j/i)j=i,..., n an intermediate partition of [a,b], 
that is, 

a = xq < x\ < ■ ■ ■ < x n -i < x„ = b, Xi-i < i/i < Xi for all i £ {1, . . . ,n}, 

while letting mesh(K) := sup i=1 n \xi — Xi-\\ go to zero. Using the Banach-Steinhaus 
theorem, one can prove that if for a right-continuous h and all continuous / the Riemami 
Stieltjes sums of (A.l) converge, h is already of bounded variation. However, we can 
weaken this assumption on the integrator by restricting the space of possible integrands. 
Recall the definitions in (2.4) and (2.5). Exploiting the concept of p- variation we now 
state an existence theorem for Ricmann-Sticltjes integrals proven by Young [17]. 

Theorem A.l. Let [a, b] C R be a compact interval, f £ 2U™ n ([a, b\) and h e 2U£ on ([a, b]) 
for somep,q > with p~ 1 + cf~ 1 > 1 . Then J f s dh s exists in the Riemann-Stieltjes sense. 

As in the classical Riemann-Stieltjes calculus, a chain rule can be proven; see [18], 
Theorem 3.1. 



Theorem A.2 (Chain rule). Let [a,b] be a compact interval and g £ 2Up 0n ([a, b]) 
for some p G (0, 2). Furthermore, let F € C X (R) with F' £ Lip(M). Then the Riemann- 
Stieltjes integral J (F' o g) s dg s exists and we have 

(Fog)(b)-(Fog)(a)= [ (F'og) s dg s . (A.2) 

J a 

At last we state a density formula, which we have not found in the literature; for a 
proof we refer to [5], Theorem 4.3.2. 

Theorem A. 3 (Density formula). Let [a,b] C R be a compact interval, f,h £ 
2H™ n (M]) and g £ W c ° n ([a,b}) for some q>0 and p > 1 with p- 1 + q- 1 > 1. For all 
x£ [a,b] we define <fi(x) := f a h s dg s . Then we have (ft £Wp° n ({a,b}) and 

b i-b 

f s d<p s = / f s hdg s . (A.3) 
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